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ABSTRACT 

It has recently been shown that the correct halo mass function high mass be- 
haviour may be obtained by an appropriate treatment of the all-mass-at-center 
problem. Here we show that the treatment of this problem leads to an accumu- 
lated mass fraction, F, which is not a universal function of the linear variance 
of the density contrast, but depends also on the shape of the power spectrum. 
For relevant power spectra F is a universal function of two quantities: the linear 
variance and a quantity related to the local shape (at the given mass scale) of 
the power spectrum. This implies, for a given shape of the power spectrum, a 
redshift dependence of the mass function, which is quite important for a cold 
dark matter power spectrum. Our results are in good agreement with numerical 
simulations. 

Subject headings: cosmology: theory — dark matter — galaxies:clusters:general 
— galaxies: halos 

1. Introduction 

The dependence of the halo mass function on redshift is a theoretical issue of considerable 
interest. On the one hand it provides a test for our set of analytical tools for understanding 
the disipationless formation of structures out of Gaussian initial conditions. On the other 
hand, it may be an important ingredient in filling in the details of the standard model. 
For example, it may be particulary relevant to obtain the prescise reionization history of 
the universe (Bromm & Larson (2004)), since the number density of objects with mass 
~ 10 5— 6 /i 1 Mq as a function of redshift must be known with some accuracy to this purpose. 
The question of the dependence of the mass function on redshift is a part of the larger question 
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of the dependence on the power spectrum, P{k). In fact, a change of redshift merely implies a 
change in the amplitude of P{k). Formaly, the accumulated mass fraction, F, is, in principle, 
a function of mass and a functional with respect to P(k): F(m, P(k)). However, the Press & 
Schechter formalism (1974) leads to mass fraction depending on mass only through the rms 
linear density contrast, a(m). In the excursion set formalism (Peacock & Heavens (1990), 
Bond et al. (1991)) there exist the possibility of an additional dependence on P(k), since, 
for top-hat filtering in real space, the relevant barrier penetration is not exactly a random 
walk process, and the way in which the actual process differs from a random walk at each 
scale depends on P(k). However, this dependence is negligible, and the same applies when 
ellipsoidal dynamics is considered (Sheth & Tormen (1999)). This implies that the mass 
function at any redshift may be obtained by rescaling cr(m); that is, that F is a universal 
function of a. 

However, evidence has been gathering from numerical simulation (Jenkins et al. (2001); 
Warren et al. (2005); Reed et al. (2006)) for an additional dependence of F on z: 

F(m, z) — F(a(m, z), z) ; cr(m, z) — D(z)a(m) 

where D(z) is the growth factor normalized to 1 at present. Now, since for a given a 
value, z is a function of m: 

F(m, z) = F(a(m, z), z(m, a)) = Ff(a, m) 

So, the additional dependence of the mass fraction on z reflects an additional (to that 
through a) dependence on mass, which, in turn, reflects a dependence on P(k), whose "local" 
shape changes with m. 

The object of this work is to show how this dependence of the mass fraction on z, and, in 
general, on P(k), may be determined using results from a recent work (Betancort-Rijo & 
Montero-Dorta (2006)), where the high mass behaviour of the mass function was explained. 
We show here that the redshift dependence is just another aspect of the same effect, and 
that the mass fraction is a universal function of two quantities, a, and a "local" spectral 
index that we represent by c. 

2. Framework 

We have shown in a recent work (Betancort-Rijo & Montero-Dorta (2006)) that the 
high mass behaviour of the cosmic mass function (at z — 0) may be fully accounted for by 
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simply adding to previous derivations (Sheth & Tormen (1999)) a consideration that has 
been, so far, overlooked. This consideration consists of the observation that, to assign a 
mass element to an object with mass to, it is more appropriate to impose the condition that 
it lies within a lagrangian distance Q(m) (Q(m) being its lagrangian radius) from the center 
of a colapsed object with mass to than to demand the mass element to be at the center of 
one of these objects. 

Using this criteria within the context of the spherical colapse approximation we found that 
the fraction of mass in objects with masses larger than to, F(m), satisfies the relationship 



exp ( — / , 2 ) 

F(m)P 2 (S l \mf > to) = _ Mm for S t > 5 C (1) 
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where P2(5i\m/ > to) is the probability distribution for the value, Si, of the linear density 
contrast within spheres of lagrangian radius Q(m) centered at a randomly chosen mass ele- 
ment belonging to colapsed objects with mass larger than to. S c is the colapse linear density 
contrast in the spherical model, and a (to) is the RMS linear density fluctuation on scale to. 
In the referred work we obtained P 2 using an approximation that allows an analytical deriva- 
tion. In this work we also use this approximation to most purposes, since it is rather accurate. 
However, with the main purpose of proving its accuracy explicity, we give here the exact 
expression for P 2 and briefly comment its relationship with that given in the reference. 

POO 

P2(Si\mf > to) = / G(mf,m)M(Si,mf,m)dmf (2) 

J m 

where G{mi,m) is the probability distribution for the mass, to/, of the object where a 
randomly chosen mass element belonging to objects with mass larger than to is to be found. 
From the definition of G(mf,m) and F(m) it follows inmediatly that: 

G(mi,m) = -—— — p — - (3) 

The quantity M{8i,mi,m) represents the probability distribution of the linear density 
contrast at a sphere with radius Q{m) centered at a randomly chosen mass element belonging 
to objects with mass mf. Through an straightforward extension of the arguments presented 
in the reference, we obtain for M: 
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M(5,,m/,m) = -^= / V ^ ^ tAfa (4) 



<r/(u) = <j/(u,Qf,Q) = (1 - 5(w) 2 )2(j(g) 

where S'(-u) is the dimensionless correlation between the linear density contrast on scale 
Q{mi) and on scale Q(m) at points lying at a lagrangian distance uQf from each other. 

1 / °° |4| 2 ^(g/fc)iy(gA;)^n(^g/)A;rfA; 
= Q, Q/) = ^ o^Qi) (5) 

3 

W(x) = —(sin x — x cos x ) 



x 3 



where \Sk\ is the power spectrum of the linear density contrast. 
Inserting expression (3) and (4) into expression (2) and this last expression into (3) we obtain 
an integro-differential equation for F(m), which, integrating by parts, may be transformed 
into an integral equation: 



F(m) 



M(5i,m,m)l y/m<r(m) 



sf 



+ / F(mt) ^-t — -dm/ (6) 

Jm dmi J 



However, in practice, it may be more convenient to use the integro-differential equation. 
Equation (6) (or its integro-differential counterpart) is valid for any 5/ > 5 C (it comes from 
expression (1)), so, for any 5i > 5 C the right hand side must be independent of 5i. 
In the referred work we implicity used M{8i,m,m) (that we called there P\(5i\5i(q = 0) = 
5 C ,Q)) in expression (2) instead of M(5i,mf, m). This implies that mass elements within 
objects with mass larger than m were assumed to lie within distance Q(m) from the center 
of the object. This is not as bad as assuming them all to be at the very center, but it is 
not as good as letting them be within distance Q(mi) from the center of the object, with 
mass ml, to which they actually belong. However, in the relevant high mass region, this 
assumption does not involve much inaccuracy, since most objects larger than Q(m) are just 
Q{m) plus a bit. 

Under this assumption, expression (2) may be transformed from an integral over mass into 
the integral of expression (4) (with w in the place of S c ) over all w > 5 C , which may be 
integrated analytically to give, when inserted into (1): 
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S(u) = S(u,Q,Q) ■ Q = Q(m) 

where S is given by (5). To simplify the evaluation of V(m), the following accurate 
aproximation for S may be used for some power spectra: 

S(u, Q, Q) ~ e~ cu2 ; c = c(Q) = Ln(S(l, Q, Q)) (8) 

c(Q) is a coefficient depending on Q but not on u, which may be obtained evaluating 
S(1,Q,Q) as indicated in expression (5). 

3. The redshift dependence 

As mentioned in the introduction, the accumulated mass function, F(m), derived within 
the excursion set formalism depends on mass only through a(m). This is referred to as the 
universality of the mass fractions, both the accumulated, F(cr), and the differential, f(a) 
(the derivative of F(a) with respect to Lna). That is, for a given value of a, the mass 
fractions are independent of the power spectrum, and, in particular, are independent of 
redshift, since changing the redshift merely implies a change of the amplitude of the power 
spectrum. However, the mass fraction obtained within the present framework (i.e. with 
an appropriate treatment of the all-mass-at-center problem) do not show this universality. 
From approximation (7), we see that F{m) does not depend on m only through a{m). V{m) 
depends on m through a and through S(u), whose functional form depends on m (through 
Q(m)) and, in general, on the whole functional form of the power spectrum. If the power 
spectrum takes the form: 



P{k) = A(k)k n ^f (k) 
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with n e ff, A slowly varying functions of k, then approximation (8) holds and V depends 
on m only through a(m) and c(m). Therefore, for these power spectra, the accumulated 
mass fraction is a universal function of a, c : 

giving a differential mass function: 

n (m) = - — ^-(a(m),c(m)) (10) 
m am 

c(m) is a sort of "local" spectral index which depends on the "local" form of the power 
spectra, but not on its amplitude. In a coarse approximation we have: 

c ( m ) ~ . k(m) ~ ^— y (11) 

For a CDM power spectrum with Q m = 0.3, h = 0.7 (BBKS (1986)) c(m) may be 
accurately obtained with the following fit: 



c(m) = 2.34740- 6 -1.678-10- 4 Ln(m)+4.888-10- 3 (Ln(m)) 2 -5.909-10- 2 (Ln(m)) 3 +0.367(Ln(m)) 4 

(12) 

Given expression (9) and (10) for the present cosmic mass function, it is inmediate to 
obtain the mass function at redshift z: a(m) must be replaced by its rescaled value a(m, z) 
while c(m) retain its value. Now, the mass corresponding to a given a value is smaller at 
redshift z than at present, so that, being c(m) a monotonicaly increasing function, the value 
of c corresponding at redshift z to a given value of a is smaller than at present. But, since, 
at fixed a, F is monotonically increasing with c, this implies that for a given a the value of 
F decreases with increasing redshift. Therefore, we see that, for a fixed form of the power 
spectrum, the dependence of F on c translates into a redshift dependence for the dependence 
of F on a: 



F(a,z)=F Q (a,c(m(a,z))) (13) 

where m(a, z) is obtained solving the equation D(z)a(m) = a. 
All we have said so far concerns the mass function obtained within the spherical colapse 
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model. In a forthcoming work we will show how to combine our formalism with ellipsoidal 
dynamics. However, the dependence of F on c arises from the appropriate treatment of the 
all- mass- at-center problem that we know is relevant only in the high mass region (a(m) < 1) 
where the spherical colapse model is a good approximation. Therefore, if we combine the 
c dependence of F obtained within the spherical model with the a dependence at z = 
obtained from numerical simulations, we may obtained a good approximation to the redshift 
dependence of the mass function. 

Let us see how this can be done. The excursion set formalism when combined with ellipsoidal 
dynamics leads to Sheth & Tormen (Sheth & Tormen (1999)) analytic mass function (with 
a = 1). In Betancort-Rijo & Montero-Dorta (2006) we conjectured that an appropriate 
treatment of the all-mass-at-center problem within the context of ellipsoidal dynamics will 
increase the S-T mass fraction by approximately the same factor as that enhancing the 
normalized P-S mass fraction when the treatment of that problem is implemented within 
the spherical model. Expression (7) tells us that this factor is simply V -1 . So, within this 
approximation we have for the present mass fraction: 



F(m) = F s _ T (a(m))(v(a(m),c(m)j) * ; F s _ T (a) = I f s _ T (af)dLna! (14) 

V J J-infty 

where F$-t stands for the S & T theoretical prediction for the mass fraction (i.e. that 
corresponding to their fitting formula with a = 1) and fs-T for its corresponding differential 
mass fraction. We have seen (Betancort & Montero, 2006) that F(m) is in good agreement 
with the fits to numerical simulations of the present mass function. 

Expression (14) implies that F(m) does not depend on m only through a. However, at 
present, a is a unique function of m, so we may find a function of a, F, such that F(a(m)) = 
F{m). Using (14) we have for F: 

F(a) = F 5 _ T (a)(v>,c(m(a))))^ (15) 

where m(cr) is the inverse function of a(m) (at present). Now, if, working under the 
wrong assumption that the mass fraction is only a function of cr, we rescaled a with red- 
shift, we would not obtain the correct mass function at that redshift, but something whose 
relationship with the actual mass function, which follows inmediately from (14) and (15), is 
given by: 



(16) 
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a(m,z) = D(z)a(m) 

By simply rescaling F(a) we are implicitly evaluating V not at the nominal mass, m, 
but at m (with cr(jn) = aim, 0) = cr(m, z)). So, to obtain the correct mass function, the 
rescaled F must be multiplied by the ratio between the values of V~ x at m and at m. For 
the mass functions themselves we have from (16) and (10): 

Mm, .) = -|7M™, «)) [ ^"'^''" gK *) + *))^^] (») 



cr(m) = D(z)(r(m) 

where f\a) is the differential mass fraction corresponding to F. Any fit to numerical 
simulations may be used for f(cr). We shall use S & T fit (a = 0.707). 
In figure (1) we show the ratio between the mass function given by (17) and the rescaled S 
& T fit to the mass function (a = 0.707). For z — 10 the ratio is close to one except for 
the largest masses. For z = 20 the actual mass function is a 35% smaller than the rescaled 
one at m — lO 5 /^ 1 M , the difference being a factor 4 for m = 10 9 /i -1 M & . For z = 30 
the difference is already a factor 2.7 for m = 10 5 /i 1 M and is larger than a factor 6 for 
m = 10 8 /i -1 M . This figure is in good agreement with the numerical results presented by 
Reed et al (2006) in their figure A3. 

Finally, to check that our results are not strongly dependent on the fact that we have used 
approximation (9), we recalculate the factor g(m,z) using expression (6), which is exact 
within the spherical model. We solve the equation iteratively using for F(mf) in the integrand 
approximation (7). We have also used for F(mf) that corresponding to the normalized P-S 
mass function (i.e. expression (7) with V — 1), obtaining very similar results. From this we 
may infer that the values of n that we have obtained with expression (6), rii, has at most a 2 
% error for all masses and redshift considered. The ratio between n\ and the value of n given 
by expression (10), n , is shown in figure (2). We see that this ratio is larger than one, thereby 
increasing the redshift dependence of the mass function, but the logarithmic correction is 
small (i.e. Ln(n )/Ln(^) 3> 1). We may therefore conclude that the results presented here 
for the redshift dependence of the mass function will remain basicaly unchanged in more 
rigorous analysis of the problem. 
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Fig. 1. — The mass function given by expression (17) normalized to the rescaled S & T mass 
function (a = 0.707) at three different redshifts. 
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Fig. 2. — Ratio between the mass functions obtained with the exact expression (expression 
(6)), m, and that obtained with approximation (17), no, at three different redshifts. 



